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Abstract 

The magnetic response of a proximity-coupled superconductor-normal metal 
sandwich is studied within the framework of the quasiclassical theory. The 
magnetization is evaluated for finite values of the applied magnetic field (linear 
and nonlinear response) at arbitrary temperatures and is used to fit recent 
experimental low-temperature data. The hysteretic behavior predicted from 
a Ginzburg-Landau approach and observed in experiments is obtained within 
the quasiclassical theory and shown to exist also outside the Ginzburg-Landau 
region. 
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I. INTRODUCTION 



A superconductor in electric contact with a normal metal induces superconducting corre- 
lations on the normal side. This proximity effect has been studied extensively, both theoreti- 
cally and experimentally (see Ref. |l| and references therein). The superconducting properties 
of the normal metal show up, e.g., in the conductance or the magnetization (Meissner effect). 
The vanishing of the resistance of a normal wire in contact to superconducting islands has 
been observed. In recent experiments (typically with wires consisting of a superconducting, 
core covered by a normal metal) by Oda and Naganoi, Mota et a/.io^ and Bergmann et alB 
the magnetization was measured. In each of these experiments, the magnetic field was ex- 
pelled from the normal-metal part (N) of the sample by superconducting screening currents 
induced by the presence of the superconductor (S). 

In the present work, we study theoretically the proximity and Meissner effect in an NS 
sandwich using the quasiclassical Eilenberger formalismiS. In materials with a high concen- 
tration of impurities (dirty limit), the quasiclassical theory can be reformulated leading to 
the Usadel equation!. In both approaches the induced pair amplitude in the normal metal, 
the density of states, the critical current, etc., have been calculated (see, e.g., Refs. |9HTTD. In 
general, for a realistic geometry the solution of the Eilenberger or Usadel equation together 
with the Maxwell equations can be performed only numerically. We have solved the com- 
bined system of equations in a wide range of parameters (temperature, external magnetic 
field, layer thickness). In addition we have paid special attention to the non-linear response 
that shows interesting hysteretic behavior at low temperatures. 

In Section II we introduce the quasiclassical Green's function formulation for the clean 
and dirty limit. The geometry and the characteristic length scales are also defined there. In 
Section III we evaluate the Meissner screening current and present some results on the space 
dependence of the various quantities. Finally, in Section IV results for the susceptibility in 
a wide range of temperatures and values of the magnetic field are presented and discussed. 
Throughout the paper we use units with Ti = ks = c = 1. 



II. THE MODEL: CLEAN AND DIRTY CASE 

Our theoretical description is based on the quasiclassical Green's function techniqueHfl. 
We will study two limiting cases, the clean limit with complete absence of impurities and 
the dirty limit with high concentration of impurities such that the elastic scattering rate 
1/Tf.i is large compared to the superconducting order parameter A and the temperature T. 
We introduce two different coherence lengths for both the superconducting and the normal 
side of the double-layer structure. In a superconducting material we have in the clean limit 



and in the dirty case 
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Here vp is the Fermi velocity, A the pair potential and D = ^vphi = the diffusion 

constant. These lengths are temperature-independent for T <C Tc. In a normal metal we 
define the temperature-dependent coherence lengths in the clean limit 



27rT 

and in the dirty limit 
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For definiteness we first consider a one-dimensional geometry with a bulk superconduc- 
tor in the region x < and a normal metal layer in the region < x < d in perfect 
electric contact with the superconductor. (Later we will also present some results for the 
cylinder geometry of the experiments.) Since we assume d ^ we can neglect the spa- 
tial dependence of the pair potential in the superconductor. This has also been confirmed 
numericallyEl. Therefore, we assume the following form for the pair potential 

A(x) = A e(-x), ImA = . (1) 

For the vector potential we use the Coulomb gauge. Together with the boundary conditions 
it is written as 



dA 

A=(0,A(a;),0), A(0) = 0, — 



= H, (2) 

x=d 



where H is the external applied field. We have neglected the magnetic field on the super- 
conducting side, i.e., we assume that the penetration depth on the superconducting side is 
much smaller than the other relevant length scales. 
We first consider the clean limit defined by 

eef(T)«/e^. 

The Eilenberger equations for this system read 



d 

{u + zeVyA{x))f3 - Afi + V:^— , gu;{v^, Vy, x) 



(3) 



which is to be combined with the normalization condition for the Green's functions 

9l{vx,Vy,x) = 1 . (4) 

The Matsubara frequencies are u = 7iT{2n + 1), fj are the Pauli matrices and Vx,Vy are 
the X- and y-components of the Fermi velocity. We assume ideal transmission, i.e. there 
are no surface potentials at the N-S-boundary and the Fermi velocities are equal on both 
sides. The boundary condition reads g^jivx^ Vy, 0—) = gujiv^, Vy, 0+). At the interface to the 
vacuum we assume specular reflection, gu}{vx,Vy,d) = g^{—Vx,Vy,d). 
The current density in the |/-direction is given by 
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j[x) = -ie2TTN{0)TY, < VyTTf3g^{v^,Vy,x) > , (5) 

where N{0) is the density of states at the Fermi level and < ■ • ■ > denotes averaging over 
the Fermi surfaces that are assumed to be spherical. 

In systems with high concentration of nonmagnetic impurities, such that 

the Eilenberger equation reduces to the Usadel equation for the isotropic part of the Green's 
function, g^{x) = < guj{vx,Vy,x) >. For our geometry it can be written as 

D-^gM-i-gM = (6) 

uj% - A(x)fi - De^A{xff2.g^{x)f2, , gui{x) 

We take the same boundary conditions as in the clean case, i. e. we neglect the magnetic 
field on the superconducting side. In the dirty limit the current density is given by the 
London-like expression 

j{x) = 'Ke'N{^)DTY.i:ihg.{x) [h , U^)] A{x) , (7) 
which together with the Maxwell equation 

= 4.,(.) = ^A(.) , (8) 
defines a local penetration depth \{x). 



III. MEISSNER EFFECT 
A. Clean Limit 

In the clean limit an analytic solution of Eq. (§) for the normal metal layer has been 
found by ZaikinS. It turns out that the fs component of the Green's function in the normal 
metal layer is spatially constant. Consequently the current density in the normal metal 
layer is spatially constant as well and can be expressed as 

Kc T /"f /"f 

1 I sin^ cos if sin 20 
~'r. Jo Jo 



47ree(Tc)3T, 



1 + (-T-)^sinh$ + ^cosh$] + cos^ 
Here we have introduced a dimensionless constant 



(9) 



K^ = 32e'N{0)vU^{T^r = -(^-^] , (10) 
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where A^v = (47re^A^(0)?T,e/?7i)~^/^ is the normal-metal penetration depth, defined analo- 
gously to the London penetration depth with the normal electron density rie replacing the 
superfiuid density. Furthermore, 

^ = —„ (11) 

vp cos u 

is the length of a classical trajectory divided by the thermal coherence length ^^(T) and 

rd 

= 2e tan ^ cosy? / A{x)dx (12) 
Jo 

is the Aharonov-Bohm phase connected with this trajectory. The last equation shows, that 
the relationship between current and vector potential is completely nonlocal. This phase 
factor leads to a shift in the energies of the Andreev levels in normal metal layeri. The fact 
that these bound states are extended over the whole thickness of the N-layer and that the 
density of Andreev levels is spatially constant, leads to a constant current density. We can 
add, that in cylindrical geometries if the normal layer is not thin comparedto the radius, 
this second condition is not satisfied. Hence, as was pointed out by Nazarovta, the current 
density is not constant in space. 

To evaluate the Meissner effect, we have to solve the Maxwell equations with the bound- 
ary condition given in Eq. (0). The solution in the N-layer is 

B{x) = H -47rjn{d-x) . (13) 

For the susceptibility of the N-layer we find 

X = . (14) 

Since the solution of the problem may lead to 47rj„c? > H the magnetic field in the region 
< X < d — H/ATT]n may change its sign relative to the applied field. This overscreening 
effect was first found by ZaikinM. In Fig. |l| we have plotted the current density and the 
magnetic field in the N-layer for different temperatures in the limit of small magnetic fields. 
Below T ~ O.lTc we find the overscreening. Below this temperature the screening quickly 
reaches its saturation value (solid curve) and a susceptibility of 3/4 as compared to complete 
screening. Thus, even at T = the screening is incomplete. 

Next we study the nonlinear response. For this purpose we write the self-consistency 
equation for the integrated dimensionless vector potential a = e A{x)dx 

1 Air 

a = e-H(f d^ejn{a) . (15) 

2 3 

Solving Eq. (|T3|) for the applied field H we get, together with Eq. (|I^), the magnetization 
curve x{H) 

TT, \ 2a ^ . . .2d 
H{a) = _ + 47rj„(a)-^ 

parameterized by the integrated vector potential. 
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B. Dirty Limit 



Here Eq. (|^) cannot be solved analytically anymore. To proceed numerically we reduce 
it to the scalar equation 



2ujF^{x) - 2AG^ - ADe'A\x)F^{x)G^{x) 



D 



G^{x)—F4x) - F^{x)—G^{x] 



(17) 



where F^^^x) = |Trfi^^(x) is the off-diagonal part of the Green's function and Guj{x) = 
(1 — Ft^(x)^)^/^ is fixed by the normalization @. The Maxwell equation is given by Eq. (H) 
with the boundary conditions given in (Q). We write for the local penetration depth 

A^^en^^S^'^""^' ^^^^ 

where we have defined an dimensionless constant = 12(Ici/'k\nY- 

We have solved this system of equations numerically for a range of parameters. Some 
results are shown in Fig. ^ for a layer thickness of 50^^ (Tc) and Kct = 100 in a weak magnetic 
field, where we neglect the term ~ in the Usadel equation. At the highest temperature, 
T = 0.3Tc, the field penetrates through the whole N-layer. Near the N-S-boundary, where the 
induced superconductivity is strong, the field decays rapidly to zero. At lower temperatures 
the field-free region increases, because the induced superconductivity extends to a region of 
linear size ~ (T). At the lowest temperature shown here, T = O.OOlTc, the field expulsion 
is London like and the field as well as the current density fall off exponentially. The reason 
is that at this temperature the local penetration depth near the N- Vacuum-boundary is 
practically constant. In the intermediate temperature regime the current fiows in a well 
defined region inside the N-layer, and consequently the screening takes place in this region. 

To calculate the nonlinear response we have to solve the complete Usadel equation Eq. (H) 
together with the Maxwell equations including the local penetration depth ([T8| ) numerically. 
From the solution of the Maxwell equations, we find the susceptibility of the normal-metal 
part 

(19) 

Results for the clean and the dirty case are discussed in the next section. 



IV. SUSCEPTIBILITY 
A. Linear Response 

First we study the case of linear response. In the clean limit we linearize Eq. with 
respect to the Aharonov-Bohm phase 0. The resulting expression has been analyzed ana- 
lytically by Zaikin in some limit^Sl and numerically by Higashitani and Nagait3. We have 
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performed the remaining integration over 9 numerically for all temperatures. In Fig. ^ some 
results for different layer thicknesses are shown. The temperature at which screening sets 
in, as well as the slope of the temperature dependence, changes drastically when the layer 
thickness increases. At T = the integrals can be solved analytically. For d ^ Xm. the 
saturation value is equal to 3/4 of a perfect diamagnet, independent of the thicknessE^. 

Here we should mention that the temperature dependence of the susceptibility cannot 
be fitted to a power law x ^ in any temperature range, in contrast to earlier theoretical 
predictions^ which suggest powers in the range 1/2 < a < 1 depending on the impurity 
concentration. Rather it satisfies an exponential law. This explains qualitatively some 
experimental resultsHi, where exponents a with values up to 2 had been found. 

In the dirty limit, we neglect the term ~ in Eq. (P) and solve the resulting differential 
equation numerically. In Fig. ^ some results of our calculation are shown. In comparison 
with the clean case, the screening sets in at higher temperatures and the saturation values 
can be larger (depending on Kd). The temperature dependence is well described by 

^ ^ [T-y^ - const) (20) 

in the intermediate temperature regime well below and above the saturation temperature. 
Our results agree with previous works in the applicable limits. In earlier wor kS , a generalized 
Ginzburg-Landau approach was used. Narikiyo and Fukuyamalll linearized Eq. (^ with 
respect to F and solved the system of equations for an infinite system numerically. However, 
our calculation is free of the limitations of the Ginzburg-Landau theory and is valid at 
any temperature. Furthermore, we have taken into account the nonlinearity of the Usadel 
equation and finite-size effects. 

The behavior shown in Fig. ^ has also been observed experimentally in very dirty 
samplesiJi. In these experiments, a temperature dependence described by x ~ (T~^/^ — 
const) was found and saturation values of the susceptibility of 90 — 95% of a perfect dia- 
magnet. 

In Fig. ^we show a comparison between theory and experiment. The experimental data0 
are measured on a cylindrical geometry with a superconducting core surrounded by a normal 
metal layer of thickness d. The theory has been generalized to this kind of geometry and also 
been solved numerically. The parameters for the theory were calculated from experimental 
values for d and lei and from theoretical values for vp and X^. No fitting parameter was used. 
The agreement is quite remarkable at low temperatures. At higher temperatures (T ~ O.IT^) 
the dirty limit theory cannot be applied anymore and this explains the disagreement in this 
temperature regime. 



B. Nonlinear Response 

We will now turn to the nonlinear response. In Fig. ^ some numerical results for the 
clean case are shown. At the highest temperature shown here, the susceptibility drops to 
zero around h ^ 0.02. At lower temperatures the solution is not unique anymore. There is 
one solution with the maximal field expulsion, the value of which depends on temperature 
(lower thick lines), and one solution, for which the field completely penetrates in the N- 
layer and there is no screening (upper thick lines). The third solution with the negative 
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derivative is unstable (thin lines). The region in which the solution is non-unique grows, as 
the temperature is lowered, but the lower boundary of this region stays at the same value of 
the field. Since only one solution can be stable, there must be a jump in the susceptibility at 
a certain point. In principle we could determine this point by comparing the free energies of 
the N-layer for the two different solutions. Unfortunately the free-energy functional proposed 
by Eilenbergeri cannot be applied for A = (normal region). 

In experiments the situation is often quite different. On lowering or raising the field, the 
jump in the susceptibility would not occur at the value predicted theoretically, but effects 
like "superheating" or "supercooling" are likely to occur. These effects were observed in 
very clean samples!. When the field is increased, the susceptibility will remain on the lower 
branch up to a certain field value and then will jump to the upper branch. On the other 
hand, when the field is decreased, the system will stay on the upper branch with complete 
field penetration. At a certain value of the field, the field will suddenly be expelled and the 
system will jump to the lower branch. These jump fields do not necessarily coincide with 
the boundaries of the instability zone. 

The boundaries of the hysteretic regions are shown in Fig. |^ for different layer thicknesses. 
The superheating field depends exponentially on temperature, whereas the supercooling 
field depends not very strongly on temperature. At low temperatures the superheating field 
tends to saturate at a certain value, which is independent of the layer thickness. Further 
calculations show that the saturation value is proportional to the material constant K^.. 

In the dirty limit there is also a temperature-dependent hysteretic effect, as shown in 
Fig. |. Above temperatures of about O.OSTc the magnetization curve is unique. Below these 
temperatures there is a region with constant lower boundary and increasing upper boundary 
as in the clean case. The absolute values of the susceptibility are drastically different. The 
saturation value at low fields and low temperatures can reach complete field expulsion for 
appropriately chosen values of Kd- On the other hand the jump in the susceptibility is less 
than in the clean case. On increasing the field the susceptibility is reduced only by ~ 50% 
at the jump and reaches slowly zero, if the field is increased further. On lowering the field, 
the jump is a slightly smaller. This hysteretic behavior was observedi. 

In Fig. y the limits of the region with hysteretic behavior for different layer thicknesses in 
the dirty limit are shown. As in the clean case the superheating field depends exponentially 
on the temperature, but not as strongly. Also, the supercooling field tends to become 
a constant at low temperatures and does not very strongly depend on temperature. In 
contrast to the clean case the saturation values of both fields at low temperatures depend 
strongly on the layer thickness. 

Finally, we would like to comment on some of our assumptions. We have assumed ideal 
interfaces between normal metal and superconductor, specular reflection at the interface 
between normal metal and vacuum, and spherical Fermi surfaces in both materials. The 
formalism presented in this paper can be modified to describe more general physical situa- 
tions. Including e.g. non-ideal interfaces will weaken the proximity effect and diminish the 
diamagnetic response of the normal metal. In this paper, we have concentrated on presenting 
a model calculation in which the main aspects of the diamagnetic response of NS sandwiches 
can be studied. Despite our simplifying assumptions, the calculation was shown to be in 
reasonable agreement with experiment, see Fig. |^. Further calculations will be necessary to 
describe the high-temperature behaviour correctly and to study arbitrary concentrations of 
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impurities, i.e., cases that are neither clean nor dirty. 

In conclusion, we have applied the quasiclassical (Eilenberger or Usadel) theory to cal- 
culate the Meissner effect in a proximity sandwich. We have evaluated the full non-linear 
magnetic response at all temperatures and for various thicknesses of the normal layer. We 
have shown that both in the clean and in the dirty limit a hysteretic behavior of the magne- 
tization of a normal metal layer in proximity with a superconductor is possible, as has been 
observed experimentally in both limits. 
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FIGURES 



FIG. 1. Current density and magnetic field in the N-layer (clean limit). Layer thickness 
d = 20^^ (Tc), dimensionless material constant Kc = 10. 

FIG. 2. Current density and magnetic field in the N-layer (dirty limit). Layer thickness 
d = 50^^ (Tc) and dimensionless material constant = 100. 

FIG. 3. Linear susceptibility (clean case) for Kc = 10. 

FIG. 4. Linear susceptibility (dirty limit) for = 100. 

FIG. 5. Comparison between experiment (AgNb)Ei and dirty limit theory. Experimental pa- 
rameters: Radius of the S-core= 35;um, d = 14.5/xm, and l^i = 2.0/xm. In the theory we have used 
d = 41.5^^ (Tc) and = 10000. 

FIG. 6. Nonlinear susceptibility (clean case) of an N-layer of thickness 20^^ (Tc) and Kc = 10 
as a function of the applied field for various temperatures. The dimensionless field is defined by 
h = H^i^{Tcf/^o. 

FIG. 7. Limits of the instability region: supercooling field h^c (thin lines) and superheating 
field hgh (thick lines) for Kc = 10. h is defined as in the previous figure. 

FIG. 8. Nonlinear susceptibility (dirty case) of an N-layer of thickness 10^^ (Tc) and K^ = 100 
as a function of the applied field for various temperatures. The dimensionless field is defined by 
h = H7TC^{Tc)y<^o. 

FIG. 9. Limits of the instability region (dirty case): supercooling field hgc (open symbols) and 
superheating field h^h (filled symbols) for K^ = 10. h is defined as in the previous figure. 
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